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We study the zero sound and the first sound in a dilute and ultracold disk-shaped normal Fermi 
gas with a strong harmonic confinement along the axial direction and uniform in the two planar 
directions. Working at zero temperature we calculate the chemical potential p of the fermionic fluid 
as a function of the uniform planar density p and find that p changes its slope in correspondence to 
the filling of harmonic axial modes (shell effects) . Within the linear response theory, and under the 
random phase approximation, we calculate the velocity c° of the zero sound. We find that also c° 
changes its slope in correspondence of the filling of the harmonic axial modes and that this effect 
depends on the Fermi- Fermi scattering length ap. In the collisional regime, we calculate the velocity 
c s of first sound showing that c s displays jumps at critical densities fixed by the scattering length 
of- Finally, we discuss the experimental achievability of these zero sound and first sound waves 
with ultracold alkali-metal atoms. 

PACS numbers: 03.75.Ss,03.75.Hh 



I. INTRODUCTION 

Effects of quantum statistics are observable with both 
bosonic [J Q and fermionic [1, 0, [E IE vapors of alkali- 
metal atoms at ultra-low temperatures. In these quan- 
tum degenerate gases, the role of dimensionality has been 
experimentally studied with bosons [E [E [l(| , while for 
fermions, there are only theoretical predictions. In par- 
ticular, it has been suggested that a reduced dimensional- 
ity strongly modifies density profiles [ill, EE EE El , col- 
lective modes [IE EH and stability of mixtures [l7l . [l8| . 
Also sound velocities have been theoretically investigated 
in reduced dimensions for both normal [IE EE] and su- 
perfluid Fermi gases [lE US Hi- 
lt has been shown that, in the weak coupling limit, 
the zero sound velocity c° and the first sound velocity 
c s of a normal Fermi gas are equal in the strictly one- 
dimensional (ID) geometry [16j . contrary to the 3D case 
where they differ by a factor V3 if the system is uniform 
[IE HE HH, or by a factor y/E if the system is cigar- 
shaped due to a isotropic harmonic confinement in two 
directions [I||. Very recently we have studied the first 
sound velocity c s of a normal cigar-shaped Fermi gas in 
the ID to 3D crossover, showing that c s exhibits jumps 
in correspondence of the filling of planar harmonic shells 

In this paper we investigate, at zero temperature, the 
zero sound velocity c° and the first sound velocity c s of 
a normal Fermi gas in the dimensional crossover from 
2D to 3D by considering a harmonic confinement in the 
axial direction. We show that by increasing the planar 
density p of fermions (or equivalently the chemical po- 
tential p) one induces the dimensional crossover in the 
system. In general, collective sound modes in a normal 
Fermi system are of two kinds. When the collision time 
r between quasiparticles is much greater than the period 
(w _1 ) of the mode itself, i.e. if lot ^> 1, the system is in 



the so-called collisionless regime. In presence of a repul- 
sive interaction, the resulting collective mode is known 
as zero sound [IE 0, HE [IE HE HI- In the present 
work we approach the problem calculating the density- 
density response function in the random-phase approxi- 
mation. The poles of this response function lead to the 
zero sound dispersion law and to its velocity of propaga- 
tion c° in the radial plane, where the gas is not confined. 
We investigate the behavior of the zero sound velocity c° s 
as a function of the chemical potential, or, equivalently, 
of the fermionic density in the transverse radial plane. 

When collisions between atoms are characterized by a 
sufficiently small collision time such that lot <C 1, the 
system may be assumed to be at thermal equilibrium. 
In this regime, the Fermi fluid is accurately described by 
quantum hydrodynamics equations where the fermionic 
nature appears in the bulk equation of state [IE HE H3 • 
In this case, the collective mode is called first sound (or 
ordinary sound, or collisional sound). The first sound ve- 
locity c s has been studied previously both in the 1D-3D 
and 2D-3D crossovers [IE HH, and in a 3D cylindrical 
geometry [32j for an ultracold Bose gas. It is important 
to stress that the propagation of ordinary sound in spin 
unpolarized Fermi gases may be studied also from the 
experimental point of view. For instance, Joseph and co- 
workers (33[ were able to excite first sound waves in an 
optically trapped Fermi degenerate gas of spin-up and 
spin-down atoms by using magnetically tunable interac- 
tions. In this paper we predict that, as in the 1D-3D 
crossover [2(| , the first sound velocity c s displays observ- 
able shell effects in the 2D-3D crossover of a disk-shaped 
ultracold Fermi gas. By increasing the planar density p, 
that is by inducing the crossover from two-dimensions to 
three-dimensions, c s shows jumps in correspondence of 
the filling of harmonic modes. Such jumps of the sound 
velocity, even if less pronounced than the ones predicted 
in the 1D-3D crossover [26|, can be experimentally ob- 
served. 
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In the last part of the paper we establish the experi- 
mental conditions to be achieved for detecting zero sound 
and first sound. In particular, by considering 40 K atoms, 
we estimate the critical mode frequency lu c which dis- 
criminates between the collisional and the collisionless 
regimes. 



II. CONFINED INTERACTING FERMI GAS 

We consider an ultracold normal Fermi gas with two- 
spin components confined by a harmonic trap in the ax- 
ial direction and investigate the regime where the tem- 
perature T of the gas is well below the Fermi temper- 
ature Tp so that we may use a zero-temperature ap- 
proach and treat separately the collisionless and the col- 
lisional regimes. The Hamiltonian of the dilute Fermi gas 
trapped by the external harmonic potential 



U(r) = ^muj 2 z z 2 , 



(1) 



and free to move in the [x, y) plane in a square box of 
size L, is given by 

H=Y. j ^( r) (-|lv 2 + ^(r)^ CT (r) 



9 -f J2 /rfr*t( r) *t /(r) ^, (r) ^ (r)(1 _ (5 _ /) 



where a = {t, |} denotes the spin variable, and the inter- 
atomic coupling is 



gp = 47T?! 5 'ap jra , 



(2) 



with ap > being the s-wave scattering length. *&o-(r) 
and ^(r) are the usual operators annihilating or cre- 
ating a fermion of spin a at the position r, therefore 
obeying anticommutation relations [H, [24|, [25|, Hz| ■ The 
field operator can be written as 



r Co 



(3) 



in terms of the single-particle operators c CTiQ (cj. a ) de- 
stroying (creating) a fermion of spin a with the single- 
particle wave function <p a (r). Because of the symmetry 
of the problem, a (r) can be factorized into the prod- 
uct of eigenfunctions of the harmonic oscillator in the z 
direction and plane waves describing free motion in the 
(x, y) plane. The field operator, then, may be written as 



**(r) 



-y 



ipj(z) exp(iK • RJCf^K 



(4) 



with R = (x,y), and where ipm(z) is the real function 



^JJ^^M-^al) . (5) 



Hj( — ) is thej-th Hermite polynomial of argument z/a z , 
a z 

a z = y/h/(mu> z ) being the characteristic length of the 
axial harmonic potential. 

The two dimensional wave vector K = (K x ,K y ) em- 
beds the translational invariance in the radial plane. By 
imposing periodic boundary conditions in the box of 
length L along x and y directions, the components K x 

and K y of K are quantized according to: K x = -y-ix 
2tt. 

and K y = ~~j~iyi where i x and i y are integer quantum 

numbers. The index a reads (j, K x , K y ), or in a more 
compact form (j, K). 

By inserting Eq. (H|) into ^ the Hamiltonian takes 
the form 



where 



H — H + H int 



a,j,K 



(6) 
(7) 



and 



Hint — 



9F 

2a 7 L 2 



y y y 



<T,<r' 3i ,h ,33 ,34 Ki,K 2 ,Q 
2 o-,ii,Ki+Q^<7' J2,K 2 -Q^ <J ' J3,K 2 C(Tj4,Ki (1 — &cj,<7') > 



with 



and 



h 2 K 2 
2m 



+ (j + \)^z , 



(8) 



(9) 



W. 



31,32,33,34 



3i,i=l V J 



hi (OH*M)H h (OHjAO- (10) 



Eq. ([9]) gives the single-particle energies of noninter- 
acting fermions. To obtain the single-particle energies 
of interacting fermions we calculate the two-spin Green 
function 

G£l K (* - = -|^(^ 1 ,k(*)4, J2 ,k(^))) (ii) 

where T denotes the time ordering operator. Here and 
in the following, we will represent the operators in the 
interaction picture (with Ho the unperturbed Hamilto- 
nian), and will evaluate the averages with respect to 
the ground state of the Hamiltonian (J6j) . In the cal- 
culation of the Green function we use a mean field ap- 
proach and replace the four operators terms involved 

in Him with (4,j 1 ,K 1 +Q £ ^4,K 1 )4' ) j 2 ,K 2 -Q^',J3,K 2 + 

c^x.Kx+Q^.Ki (4^ 2 ,K 2 -cA',i3,K 2 ). The poles of the 
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FIG. 1: (Color online). Chemical potential p vs. planar 
density p. Three values of the scaled interaction strength: g — 
10~° (dotted curve), g = 0.35 (dashed curve), g = 0.6 (solid 
curve). The chemical potential is in units of hto z , density in 
units of and lengths in units of a z . 



Green function with a = a' and ji = j2 = j give di- 
rectly the single-particle energies £j.k of the interacting 
fermions 



€j,K = -7T— + U + o) huJ z + 



2m 



2a z L 2 



E 

1=0 



Uj,ffll 



where 



K.a 



(12) 
(13) 



is the total number of fermions in the harmonic shell I 
and Uj t i — ui t j — Wi t i : j t j is obtained from Eq. (fT0|) . 

In Tab. 1 we report the numerical values of the coeffi- 
cients uij for I and j up to 4. 

The maximum positive integer l max involved in the 
summation at the right hand side of Eq. (fl"2")l clearly 
depends on the single-particle state (j, K) and on the 
interaction strength gp. 

If spin-up and spin-down states are equally populated 
we may write the total number of fermions in the i-th 
harmonic state as 



(14) 



K 



with p, the chemical potential, i.e. the Fermi energy. 

From here we may calculate the 2D fermion density in 
the (x,y) plane as a function of the chemical potential. 
Under the condition L a z on the planar box size, the 
summation over K may be repaced by an integral to get 
the 2D density associated to the jth axial mode, which 
is uniform in the (x,y) plane and is given by 



L 2 (2tt) 2 



d 2 KQ(fl-e^ K ) . (15) 



(U) 




(0,0) 


1 


(0,1) 


1/2 


(0,2) 


3/8 


(0,3) 


5/16 


(0,4) 


35/128 


(1,1) 


3/4 


(1,2) 


7/16 


(1,3) 


11/32 


(1,4) 


75/256 


(2,2) 


41/64 


(2,3) 


51/128 


(2,4) 


329/1024 


(3,3) 


147/256 


(3,4) 


759/2048 


(4,4) 


8649/16384 



Table 1. Coefficients uij. Note that it;, 



If only only the lowest (j — 0) single-particle mode 
along the z axis is occupied, then the system may be 
considered as strictly 2D. After inserting Eq. (fT2|) into 
Eq. (fT5|) an easy integration gives 



where 



nti 



9F_ 

2a z 



2 



1=0 



(16) 



(17) 



is the chemical potential measured with respect to the 
ground state energy of the axial harmonic potential. The 
expression (|16p may be rewritten as 



E ( 5 3,l+9Uj,i)pl 



1=0 



m 



(p - huj z j) 



where 



a 



9f 



2nhuj z 



2a_F 



(18) 



(19) 



is the scaled inter-atomic strength. 

The total planar density is then given by 

'max 

p = £ pi > 



(20) 



where the only nonzero densities pa, pi, pi mam -i, pi mam 
are obtained by solving the l max + 1 Eqs. (fT%|) at fixed p, 
and g. The value of l max is determined by requiring that 
all the l max + 1 solutions are positive. Note that when 
only p is non zero, we have a strictly 2D system. In 
Fig.l we plot the chemical potential p as a function of the 
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planar density p for three values of the scaled interaction 
strength g. The figure shows the changes in the slope 
of these curves which occur when the planar density is 
sufficiently large so that a new excited level of the axial 
harmonic well (a new shell) begins to be populated. 



III. ZERO SOUND 

In the present section we analyze the collective mode 
of the fermionic gas in the collisionless regime, i.e zero 
sound. To this end we calculate the density-density re- 
sponse function: 

x a ' CT (Q,t) = 

~e(t)<[MQ,t),ftUQ,o)]> ( 21 ) 

since its poles provide the dispersion relation and the 
velocity of propagation of this mode. 

In Eq. (f2"Tj) we have introduced the 2D density number 
operator ^(Q) according to 

*V(Q) = E *V,j,K(Q) = E '•„.., |< r " .. K - (J • (22) 



At any time, we may define 

KQ,*) = I>a,i,K(Q,*) 



(23) 



K 



We proceed by calculating the equations of motion for 
spin diagonal [a — a') and non diagonal (er ^ a 1 ) density- 
density response functions 



= 8(t) < K J , K (Q,t),<,(Q,o)] > 

- ~e(t) < [PW(Q,t), #(*)], 4j(Q.°)] >, 

= ~0(i) < [[n CT)J - K (Q ) t),H(t)],nt, > .(Q,0)] > . 

After a len gthy but straightforward calculation and us- 
ing RPA [2J, l28| to decouple the higher order response 
function coming from the commutator on the r.h.s. of 
the equation of motion, the time Fourier transforms 
Xn(Q) u ) 01 the spin diagonal and non-diagonal density- 
density response functions read, respectively 

x£f(Q> w ) = xB{Q,,u) 



a z L 



2 a jj 



1=0 



1 !0) (Q>w) ^ Mj,/x^ ,,T (Q,w) 

' z 1=0 



The function xj°-(Q,w) = 
both equations, is given by 



(0),cr,cr \ 

j,j w~ J ~ XjJ (Q' w )' appearing m 



xS(Q.«)=E- 



n<r,j,K. — «<tj,K+Q 



^hu- (e° K+Q - e° K ) + iftjj 



(24) 



where the positive imaginary term in the denominator 
ensures causality and the limit 77 — > + will eventually be 
taken. The totally symmetric density-density response 
function defined as 



1 



[x T ' T (Q, w) + x i,T (Q, w) + x u (Q, w) + x u (Q, w)] 

(25) 



finally reads 

Xjj(Q> w ) 



xgj(Q^) 
2 

5F (0) 



+ TT2Xj'°i(Q»w) ^ u^xwCQ. 



J=0 



(26) 



This equation is one of the main results of the paper. 
It may be recast in the form of a system of equations 
for the response functions Xz,i(Q' w ) pertaining to the 
various axial modes I: 



E K< 

1=0 v 



,9F 



(27) 

in terms of which the full response function is given by 
X(Q,w)= Exm(Q,w) (28) 

The poles of this response function coincide with the 
zeros of the determinant of the matrix of coefficients of 
X;,i(Q,oj) in the system of equations (j2"T|) in the limit 
X] — ► + . Since zero sound is appreciably damped for 
large Q, (24|, we look for zeros of the determinant in the 
limit Q — ► 0. In this regime, one may evaluate x\°i (Q> w ) 
of Eq. (f2"4"| under the condition 



to find 



S= Q > V 2 ^ 



mL 



(29) 



(30) 



where use of Eq. (fTo) has been made. Note that 



nh 2 



VFJ 



-pi 



(31) 
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is the planar Fermi velocity in the Z-th axial shell (I — 
0, 1, l m ax)- v f,i corresponds to the Fermi velocity of 
a strictly 2D uniform Fermi gas with density pi, which 
at fixed p is & function of g. Using Eq. (|28[) . we may 
calculate the response of the system for any value of the 
chemical potential, i.e. for any number of excited levels 
of the axial harmonic trap. We start by considering the 
case when only the lowest state of the harmonic axial 
well is populated, i.e. l max = 0. Then, from Eq. (j27|) . 
one gets the density response function XofiiQ^) as 



X$(Q,«)/2 



Xo,o(Q,w) = 



1 - (9F/a z L 2 )u Q fiXol(Q,u 



(32) 



By using Eq. f30J| with I 
poles of Eq. (|3"2")l satisfy u> - 
zero sound is given by 



0, we find easily that the 
c° s Q where the velocity of 



c° = v F ,o f(g) 



with 



/(<?) = 



i + g uq,o 

y/l + 2 g u ,o 



(33) 



(34) 



the term which takes into account the effects of interac- 
tion. We remind that the chemical potential p, is related 
to the 2D density p by Eq. (fig)), i.e. 



■nh, 



(1 + .9 «o,o) P ■ 



(35) 



The expression (|3"31 for the zero sound velocity coincides 
with that of a strictly 2D Fermi gas 29] of density po, 
with the effect of the axial harmonic trap embodied in 
the definition of the dimensionless parameter g uq,q as 
the effective 2D coupling. We recognize in this effective 
2D coupling the familiar Landau parameter F$ [2^] for a 
strictly interacting 2D Fermi gas. Let us denote g ito,o 
by Fq . By using this definition and the Eq. (|35p , we may 
express the zero sound velocity ([33]) in terms of F 




FIG. 2: (Color online). Zero sound velocity c° vs. planar 
density p. Three values of the scaled interaction strength: 
g = 10" 5 (dotted curve), g = 0.35 (dashed curve), g — 0.6 
(solid curve). The zero sound velocity is in units of a z cu z , 
density in units of l/a 2 z , and lengths in units of a z . 
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FIG. 3: Ratio (Q/vf,o vs. chemical potential p. Three values 
of the scaled interaction strength: g = (dotted curve), g = 
0.35 (dashed curve), g = 0.6 (solid curve). Note that the 
dotted curve is superimposed to the solid one. The velocities 
Cg and vf,o are in units of a z uj z , lengths in units of a z , and 
the chemical potential in units of hco z . 



2 -(, 

m 



Trhuj z alpFo) 



1 + F 
VI + 2 F ' 



(36) 



Now we consider a density p large enough so that some 
excited axial harmonic levels are populated. In this case, 
the velocity of zero sound is given by the zero of the 
determinant of the matrix in the l.h.s. of Eq. (|27[) with 
lu > Q. The calculation is again straightforward and the 
results are reported in Fig. 2, where we display the zero 
sound velocity c° as a function of the planar density p. 
As it happens in the strictly 2D case, for a given value 
of the planar density p, c° s increases by increasing the 
repulsive two body interaction g. 

In Fig. 3 we plot the quantity c®/vf,o as a function 
of the chemical potential p,. When g = one finds that 
Cg=VFfi. If g > but p, is sufficiently small so that only 
the lowest axial mode is populated, then the quantity 



c®/vf,o does not change when the chemical potential is 
varied (dashed line). This behavior corresponds to the 
one described by Eq. (|33[) . When the chemical potential 
is greater than a critical value p c , excited axial modes 
begin to be occupied by the fermions and c®/vf,o grows 
with p (see the inset of Fig. 3). 



IV. FIRST SOUND 

In this section we analyze the collisional case of a nor- 
mal paramagnetic Fermi gas. To analyze a sound wave 
that travels in a planar direction it is important to deter- 
mine the collision time r of the gas. According to and 
[HI, if the local Fermi surface is not strongly deformed 
then t = t (T/T f ) 2 , where r = 1 / '(p (3D) av F ) with a 
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the scattering cross-section and p^ 3D ^ the 3D density. In- 
stead, if the local Fermi surface is strongly deformed, as 
in our case, then the collision time is simply given by 
t = To [Hj]. It is easy to find that in our 2D geometry 
tq ~ l/{uj z pa 2 F ). Setting u the frequency of oscillation of 
the sound wave, the wave propagates in the collisionlcss 
regime if lot 3> 1 and in the collisional regime if lot <C 1. 

In the collisional regime by adopting the equations of 
hydrodynamics [23|, [25|, H3] the planar first sound veloc- 
ity c s of the Fermi gas is given by the zero-temperature 
formula 



(37) 



Cs = \ / — 




By using this formula and Eq. (|18p we determine the 
behavior of c s versus p and of c s /vf,o as a function of p, 
as shown in Fig. 4 and in Fig. 5, respectively. These 
figures, obtained both for g = and for g > display 
shell effects, namely jumps of the first sound velocity c s 
when the atomic fermions occupy a new axial harmonic 
mode. 

If the interaction is negligible in the equation of state, 
the values p c of the density at which such jumps occur 
may be analytically obtained from Eq. (fl~8]) . In fact, we 
obtain the behavior of the planar density as a function 
of the chemical potential by setting g = in the Eq. 
(|18p and performing the sum over the quantum number 
j. Then, we get 



j'=o 



1 



nai 



JL 

TlLO z 



irai hu z 



where 



S(x) 



(1 + I[x])(x--I[x]) 



(38) 



(39) 



with I[x] the integer part of x, Eq. (J3SJ) with Eq. ([39]) 
represents a simple analytical formula which gives the 
uniform planar density p of the Fermi g function 
of the chemical potential p.. The values p c signing the 
jumps can be obtained from Eq. ([38]) with with Eq. ([39]) . 
The jumps of the first sound velocity are given by 

Po= fc( fc+ 2 1} with fc= 1,2,3,.... (40) 

When the critical value of the density is approached from 
the left, from Eq. (|37|) we find that the sound velocity is 



c, = LO z a z 



k + 1 



(41) 



while when n c is approached form the right, we find 
+ [k 



(42) 



The height of the first jump of c s is then given by 



Ac. = . 



V2 



(43) 



When interactions are not negligible in the equation of 
state, one must use Eq. ([T8]) to get the values of p for 
which the slope of p exhibits discontinuities, and calcu- 
late the corresponding discontinuities in c s . For example, 
at the first discontinuity, which appears when l max goes 
from to 1 the value of the density is 



1 



1 



1 1 



ira 2 1 + gu 0fi Tral 1 + F 



(44) 



The corresponding discontinuity exhibited by the first 
sound velocity is 



where 



c s = lo z q z] 



(1 + fftio,o)[l + ff(«o : o - 



(45) 



(46) 



and 



u z a z 



1 + g[u 0:0 + ui,i + .g(uo,owi,i - u 2 o,i)) 



[1 + g(u 0fi - «i,i)][2 + g(u 0fi 



2u ,i + 

(47) 



where the it's are given in Tab. 1. We observe that these 
shell effects are similar to the ones occurring in the 1D- 
3D crossover (2(|, where the jumps are sharper because 
the sound velocity goes to zero in correspondence of the 
filling of planar harmonic modes, showing the crucial role 
played by the dimensionality. In the next section we shall 
discuss the experimental conditions necessary to observe 
shell effects in the velocity of first sound. 

Now we comment on the 2D-3D crossover in our 
fermionic gas, starting from the strictly 2D case, cor- 
responding to lmax = in Eq. ([T6]) . In this case the 
chemical potential p is given by Eq. (|35|) and therefore 
from Eqs. ([3T]) and ([3"T|) one finds 



VF,0 

V2 



\J (1 +5Mo,o) 



(48) 



or, in terms of the Landau parameter Fq 



^-(p-nhLo z a 2 z pF ) 



V2 



(49) 



recovering the usual result for a 2D Fermi system. Such 
a behavior is illustrated by the first plateau in Fig. 5. 
Incidentally, note that when Fq ~ the zero sound ve- 
locity ([3"6"]l is larger than the first sound velocity (14"9"]) by 
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FIG. 4: (Color online). First sound velocity c s as a function 
of the transverse density p. Three values of the scaled interac- 
tion strength: g — (dotted curve), g — 0.35 (dashed curve), 
g = 0.6 (solid curve). The first sound velocity is in units of 
a z ui z , density in units of l/a^, and lengths in units of a z . 
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FIG. 5: (Color online). The ratio c 3 /vf,o versus chemical 
potential p. Chemical potential is in units of Twj z . We report 
behavior for g = (dotted curve), g = 0.35 (dashed curve), 
g = 0.6 (solid curve). The velocities c s and vf,o are in units 
of a z ui z , lengths in units of a z , and the chemical potential in 
units of hu z . 



a factor of \/2. In our case, the quantity Fo is smaller 
than one. Then the two velocities are different between 
each other also when Fq remains finite. 

We observe that the first sound velocity depends ex- 
plictly on the chemical potential. In fact, remembering 
the relation between p and p given by Eq. (|18p . and us- 
ing Eqs. (|38[) and (f3"9")l . one may show that the velocity 
of first sound takes the form: 



We observe that starting from the familiar Thomas- 
Fermi formula for the 3D local density p(z) of the Fermi 
gas under axial harmonic confinement, given by 



p{z) 



(2m 



,3/2 



2 2 

muj, z 



.3/2 



3tt 2 ^ 3 

and integrating over z, one obtains exactly Eq. 



(54) 




_ 1 I[p,/fruj z 
2 pi/Tiuiz 



tig) 



p/huj z 



(50) 



where t(g) is a complicated but analytical function of g 
and of the matrix elements u's, which vanishes if g = 
0. We stress that the chemical potential p depends on 
the planar density p and on the interaction strength g. 
Fig. 5 shows how c s /vp t o explictly depends on p when 
axial harmonic states other than the ground states are 
occupied. From Eq. (|50[) one finds that for any g at very 
large p (i.e. also very large p) the speed of sound is given 
by: 



(51) 



In absence of two-body interactions (g = 0) for p > 
l/(ira 2 ), i.e. for p, > Tiuj z , several single-particle states of 
the axial harmonic well are occupied and the gas exhibits 
the 2D-3D crossover. Finally, for p ^> l/(Tra 2 ), i.e. for 
p, ^> huj z , the Fermi gas becomes 3D. In this case from 
Eqs. (J37J), §B and ([39]) one finds 



p=(27r) 1 / 2 hu Jz (alp) 1 / 2 



= V2 (27T) 1 



(alp) 



1/4 



(52) 



(53) 



V. EXPERIMENTAL FEASIBILITY 

We now discuss the experimental conditions to detect 
the zero and first sound in a quasi-2D two-component 
Fermi gas. We suppose that in the [x, y) plane the system 
is in a square box of length L = 20 /im, while in the axial 
direction z there is a harmonic confinement characterized 
by a frequency of 160 kHz, which corresponds to a char- 
acteristic length a z ~ 0.1 pm. We consider 40 K atoms 
trapped in the two hyperfine states |/ = 9/2, mt = 7/2) 
and |/ = 9/2, mt = —7/2) with / the total atomic spin 
and m/ its axial projection. In this case the s-wave scat- 
tering length reads of = 174 ao, where ao = 0.53 • 10~ 10 
m is the Bohr radius (35|. For this value of ap, we can 
determine the value of the collective mode frequency u> c 
which discriminates between the collisionless and the hy- 
drodynamic regimes. In fact, the collision time r can 
be estimated as r ~ 1/ (oj z pa F ) and the corresponding 
critical frequency is uj c — 2tt/t ~ ui z (2irpa F ). If the fre- 
quency uj of the collective mode is much greater than uj c 
(i.e. for uj oj c ) the system is in the collisionless regime; 
otherwise (i.e. for uj <C uj c ) the collective excitation is 
hydrodynamic-like. 

One can experimentally change uj c in the Fermi gas 
of 40 K atoms by varying the planar density p. For in- 
stance, from Fig. 4 and Eq. (39) one finds that in corre- 
spondence to (pa 2 z )n = 1/2 the system is strictly 2D. In 
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this case, the planar density is p ~ 1.6 ■ 10 13 atoms/m 2 , 
the total number of atoms is N = pL 2 ~ 6.4 • 10 3 , and 
the critical frequency reads lo c ~ 1.4 kHz. In addition, 
the zero sound velocity is = i>f,o — 0.02 m/s, while 
the first sound velocity is c s = vf,o/V% — 0.014 m/s. 
Moreover, one finds that the scaled interaction strength 
is g = 2ap/a z ~ 0.18. To conclude, we observe that one 
can change the critical frequency lo c also by varying the 
scattering length oj? via magnetic Feshbach resonances. 

VI. CONCLUSIONS 

In the present paper we have considered a dilute and 
ultracold interacting disk-shaped Fermi gas confined in 
the axial direction by a strong harmonic trap and uni- 
form in the two planar directions. In particular, we have 
analyzed the behavior of the velocity of propagation of 
the zero and first sound modes in this gas. 

In the first part of the paper, we have studied the zero 
sound mode from the density-density response function 
within the linear response theory and under the random 
phase approximation. We have obtained the zero-sound 
dispersion law by calculating the poles of the density- 
density response function. Then we have analyzed the 
behavior of the velocity of the zero sound as a function of 
the planar density for different values of the Fermi- Fermi 
scattering length. We have verified that as in the strictly 
two dimensional case, in correspondence of a given value 
of the planar density, the zero sound velocity increases 
as a consequence of the increasing of the strength of the 
fermion-fermion interaction. We have studied, moreover, 
the behavior of the ratio between the zero sound velocity 
c~ and the planar Fermi velocity vf,q at varying of the 
chemical potential; when, in presence of a non vanish- 
ing fermion-fermion interaction, the chemical potential 
is greater than a critical value, excited axial modes begin 
to be occupied by the fermions and <P s /vf,o grows with 
fi. In principle, the same kind of analysis may be carried 
out in presence of a harmonic potential in the transverse 



radial plane instead of the axial direction (see [26j). In 
this case, the treatment of the problem becomes much 
more complicated. In fact, the cigar-shaped configura- 
tion is characterized by two quantun mumbers related 
to the harmonic energetic levels in the planar directions. 
However, we expect for the zero sound velocity the same 
behavior as in the case of the only axial harmonic trap- 
ping. 

In the second part of the paper, from the formula which 
relates the chemical potential of a disk-shaped Fermi gas 
to its uniform planar density, we have calculated the col- 
lisional sound velocity of the system. We have found 
that this sound velocity gives a clear signature of the 
dimensional crossover of the two-component Fermi gas. 
Our calculations suggest that the dimensional crossover 
induces shell effects, which can be experimentally de- 
tected. In fact, as in 1D-3D crossover (26|, also in the 
2D-3D crossover the sound velocity exhibits jumps in cor- 
respondence of the filling of axial harmonic modes. We 
have investigated the behavior of the ratio between the 
first sound velocity c s and the planar Fermi velocity as 
a function of the chemical potential. We stress that our 
study shows that when only the ground state of the har- 
monic well is populated and when the chemical potential 
is sufficiently large so many excited axial modes are oc- 
cupied, our results reproduce the behavior of a 2D and a 
3D Fermi system, respectively. 

Finally, we have discussed the possibility of observing 
in experiments the zero and first sound by using gases of 
40 K atoms trapped in the two lowest hyperfine states. In 
particular, we have estimated the sound-mode frequency 
which discriminates between the collisionless and hydro- 
dynamic regime. We have found that these regimes re- 
quire severe geometric and thcrmodynamical constraints, 
but they can be reached with the available experimental 
setups. 

This work has been partially supported by Fondazione 
CARIPARO. The authors thank Giovanni Modugno for 
useful suggestions. 
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